A1 = readmatrix('NiI2_20Vcool_30K.csv')
A2 = readmatrix('NiI2_-20Vcool_30K.csv')
B = readmatrix('supercontinuum_power_calibration.csv')

m = B(:,1)
n = B(:,2)*0.1

x1 = A1(:,1)
y1 = A1(:,3)*1e12./n
z1 = A1(:,3)*1e11/2

x2 = A2(:,1)
y2 = A2(:,3)*1e12./n
z2 = A2(:,3)*1e11/2


hold on
xlim([450, 800])
plot(x1,y1,'Linewidth',1,'Color',[0.8 0 0])
plot(x2,y2,'Linewidth',1,'Color',[0 0 0.8])
%plot(x1,z1,'Linewidth',1,'Color',[0.8 0 0])
%plot(x2,z2,'Linewidth',1,'Color',[0 0 0.8])
xlabel('Wavelength [nm]')
ylabel('Photocurrent [pA/uW]')
ylim([-0.2,0.2])
set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
legend('2 MV/m', '-2 MV/m')
legend boxoff
box on
hold off

set(gcf,'Position',[0, 0, 212, 200])